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Abstract 


In this work we studied the acoustic variability in the induced release call on wild populations of the European common 
frog Rana temporaria. We used a computational model to analyze the induced release calls in both sexes. Our analysis 
detected three different acoustic phenotypes: (a) Non-canonical males comprised a large part of the population and were 
the main cause of the biased sex ratio (Phenotype 2), since (b) canonical males (Phenotype 1) and (c) canonical females 
(Phenotype 3) had close to a 1:1 ratio. Phenotype 2 were much more abundant in locations were tadpoles were exposed 
to high water temperatures during their development, suggesting that a masculinizing effect of temperature is 
responsible for the biased sex ratio, and a likely origin for “pirate” males which search for newly laid clutches to fertilize 


not fertilized eggs by the “parental” males. 
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Introduction temporaria, such as a wider skull, a lower distance 
between eyes and a shorter and wider mouth, differences 

The European Common Frog is the most widespread in the femur sinuosity,...Later, Teacher, et al. [2] studied 
anuran species of Europe inhabiting all kinds of the distribution of the European clades of Rana 
permanent and temporary waters from northern Spain to temporaria and claimed that there were two different 
subarctic Fennoscandia, and all the way to the Ural clades in the Pyrenees which overlapped our study area 
mountains in the east [1]. in the Central Pyrenees. Taking into account that one of 


the main factors acting in the call characteristics is the 
In 1995, the third author, described a morphotype skull dimensions [3]. 
which inhabit the Central Pyrenees designed as Rana 
temporaria aragonensis which presented morphological In 2004, the third author, found behavioral differences 
differences when compared to Rana _ temporaria in the mating strategies of our study population, where a 
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previously unknown reproductive behavior in vertebrates 
with external fertilization called ‘clutch piracy’ was 
described in our study area with a male-female sex ratio 
2:1 [4]. This “pirate” males search for recently laid 
clutches clasping them to fertilize the eggs not yet 
fertilized by the “parental” male. 


It is possible that this longstanding sex imbalance has 
been exacerbated by male aggression resulting from a 
lack of available females [5,6]. We have in several 
occasions observed males attacking and kicking females. 
In addition, this sex imbalance could be the cause for 
novel mating strategies such as the clutch piracy 
previously reported in this population [4]. This kind of 
extreme environmental masculinization can be dangerous 
for populations. The combination of unbalanced 
phenotypical sex ratios, reduction in genetic variability 
and male aggression are likely to result in drastic 
population reduction or even extinction [7]. The possible 
impact of temperature rise on animals with temperature- 
based sex determination is not a new idea. It has been 
observed that even small increases of 1-2 °C in air 
temperature could have significant effects on the sex 
ratios of sea turtles, and some authors warned about 
large, fast temperature rises having very detrimental 
effects on populations [8]. Some authors Miller, et al. [9] 
argue that it is possible that dinosaurs had temperature 
sex determination and that fast climatic change at the 
Cretaceous-Tertiary boundary was a contributing factor 
to the collapse of their populations. Lastly, recent studies 
Mitchell, et al. [10] Grayson, et al. [11] point out that 
species that show female to male sex reversal at high 
temperatures such as  Sphenodon tuatara_ could 
experience severe impact under current predictions for 
temperature rise. Our data shows that the existence of 
more intense and frequent high temperature extremes 
can have a severe effect on isolated populations that show 
temperature-mediated masculinization. Populations with 
a sex ratio heavily skewed towards males such as the one 
reported in this study are at risk of demographic collapse 
[11]. 


Although masculinization at high temperatures is well- 
known in R. temporaria as well as other amphibians [12], 
it is unclear exactly how widespread it is. If this is a 
common phenomenon among amphibians, temperature- 
mediated demographic collapse could be an important 
factor in the declines observed in the past years [13,14]. 
This is consistent with the observation that declines have 
a tendency to occur at higher altitudes [15], as shorter 
activity periods would delay tadpole metamorphosis to 
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the warmer summer months, as described in our study 
population. 


Our main goal is to find the existence of acoustic 
phenotypes, in wild environments and not in laboratory, 
due to the occurrence of both clades of R. temporaria in 
our study area. Indeed, another objective is to acoustically 
differentiate “pirate” from “parental” males in our 
population. Lastly, the third objective is to study the effect 
of temperature during metamorphosis in the acoustic 
phenotypes. 


Material and Methods 


Study Area 


During the summer of 2013 and 2014 we sampled 68 
water bodies in an area of 600km? located in the Central 
Pyrenees (Figure1). From which we selected seven 
sampling sites where we found adult frogs with a snout- 
vent-length (SVL) between 4.6 and 6.4 cm (approx. 5 to 6 
years old). Three spring and streams located in low- 
altitude sampling sites, Aguas Tuertas, Formigal, 
Surgencias en Acherito (valley), and four glacial lake high- 
altitude sampling sites, Bachimafia Bajo and Alto, Brazato 
and Asnos. Both with an average distance higher than 15 
km and large geographical barriers. 


Figure 1: Map of our study area in the Central Spanish 
Pyrenees (Aragon, Spain) with the sampling sites 
marked by blue dots. White squares placed every 3 
minutes North-West. Green and red level lines 
represent altitudes at 1600 and 2400 m ais. 
respectively. Black dash line is the Spanish-French 
border. 


Animal Sampling 


We sampled a total number of 532 frogs from which 
only 55 were sexually receptive males and females, 29 in 
springs and streams and 26 in glacial lakes with altitudes 
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ranging from 1600-1700 and 2200-2300 m aiszl. 
respectively. We identified sex unambiguously using 
external morphological sexual characters (nuptial pads, 
size, and white throat in males). These frog populations 
breeded every year. 


In order to avoid the size effect in our samples, we 
selected the induced release calls of adult frogs with a SVL 
in a specific range 4.6 to 6.4 cm. The SVL data collected 
from smaller frogs were not considered for this study. 


Sound Data Recordings 


Although other European populations of R. temporaria 
have been described to sing in chorus during the breeding 
season [16], individuals from both sexes of R. temporaria 
in our alpine study population are known to do not sing in 
chorus during the breeding season and only emit 
spontaneous vocalizations when they are underwater, 
which led us to use induced, above-water release calls for 
being less complicated to obtain [3]. 


We recorded induced release calls of individuals from 
our study population and analyzed them in order to test 
for the presence of different acoustic phenotypes. Release 
calls are produced by both male and unreceptive female 
frogs when grasped by male frogs [17]. This type of 
vocalization can be easily induced by gently grabbing the 
frog and results in little to no damage or stress for the 
animal. 


Vocalizations were recorded using a TCD-D8 Portable 
Digital Audio Tape (Sony Corporation, Tokyo, Japan) a 
directional microphone (C568B, AKG acoustics, Vienna, 
Austria) with integrated amplifier and 20-20.000 Hz 
bandwidth. Sampling rate was 44.1 kHz with a sampling 
depth of 16 bits. 


Our recordings consisted of release calls, which are 
produced by both male and unreceptive female frogs 
when clasped by male frogs [18]. This type of call 
presents sexual dimorphism in the number of pulses, 
inter-pulse interval, maximum frequency and duration of 
the call [19,20]. This call can also be used to distinguish 
Neotropical and European toads’ species [21,22]. 


The release call in R. temporaria was classified by van 
Gelder, et al. [23] in function of the amplitude pattern of 
the pulses. We studied the pulses defined as a “vibrator 
note” in function of the amplitude pattern of the peaks, 
which is the main difference with respect to the other 
bioacoustics amphibian studies (Figure 2) [3]. These 
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pulses with peaks are an inferior fractal order with 
respect to the entire call. 


E 00:00:00 00:00:02 00:00:04 


00:00:00,000 


00:00:00,050 


E 00:00:00,000 00:00:00,010 


6.0 


Figure 2: A Oscillogram showing the different pulses 
of a entire call. B Oscillogram of a specific pulse 
extracted from the A, showing the different peaks in it. 
C Oscillogram of a specific peak extracted from B 
showing its amplitude pattern. These pulses with 
peaks represent an inferior fractal order with respect 


to the entire call. 


Environmental Data 


We measured the water temperature with a digital 
thermometer in the majority of the sampling sites during 
the end of July and the first week of August, under the 
same climatic conditions, sunny, and between 12 and 14 
solar hours. 


To be able to extrapolate the water temperatures 
along the year we calculated the theoretical solar 
radiation in our sampling sites using the Morgen, et al. 
[24] for the parallel 50 (close to the study area), which 
already consider the slopes and shielding effect of slopes. 
In each sampling site we traced two profiles using the 
Global Mapper GIS application (Company: DroneMapper) 
software, one North-South profile and a second West- 
East. Afterwards we calculated the profiles angles with 
the highest slope in the sampling site and its exposition 
using the GeoGebra International software [25]. 
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Data Processing 


For each of our analyzed call pulses (Figure 2) we 
characterized the following variables that define some of 
the fundamental characteristics of a vocalization using 
custom Matlab scripts (MATLAB R2013b 8.2.0.701, 
Mathwoks, INC) to repeat the analysis of 2152 pulses 
from the 55 release calls emitted by the 55 adult frogs: 


e The first harmonic or fundamental frequency (FO). 
e Number of peaks in each pulse (Figure 2). 
eA measure of distortion calculated using the thd 
MATLAB function. The total harmonic distortion is 
determined from the fundamental frequency and the first 
five harmonics using a modified periodogram of the same 
length as the input signal. It is a negative number which 
shows the frequency modulation, being lower as the 
modulation increases. 

In addition, added other variables 


we eight 


corresponding to the location, year, month and day of 


A 


Canonical discriminant function 


Function 2 


Function 1 


Figure 3A: Results of the discriminant analysis that resulted in the classification of pulses in four types. X and Y axes 
show the value of Functions 1 and 2, respectively. Cluster centroids are marked by a square and its group number. 


3B: Classification of the studied vocalizations in four classes of call according to the proportion of each type of pulse 


present. 


Ratio of pulses types 


LL 


sampling, snout-vent length (SVL) in cm, a consecutive 
number of an analyzed pulse within the call and sex of the 
individual as determined by morphological 
characteristics. 


Using the IBM SPSS Statistics 19.0 (SPSS, Inc., Chicago, 
IL) statistical package for Windows, we carried out a 
hierarchical cluster analysis (HCA) using the FO, Number 
of peaks in each pulse and Distortion variables clustering 
2152 pulses into four clusters. 


Using each cluster as a grouping variable and the 
above-mentioned variables as independent variables we 
performed a discriminant analysis, determining as new 
clusters the predicted groups for each pulse and using 
them as a base for the next discriminant analysis. This 
process was repeated three times until the clusters were 
fully stabilized, thus optimizing the discrimination, each 
of these stabilized clusters are called Type of pulse 
(Figure 3A). 


B 


BPulse type 1 
Pulse type 2 
Pulse type 3 

BiPuise type 4 


Class 1 Class2 Class3 Class4 
Classes of call 


Given the four different types of pulses, each call (one 
call for frog) was subsequently classified using a cluster 
analysis taking as variables the ratio of pulses types 
present within it. We obtained four classes of call (Figure 
3B). To measure the strength of these four variables 
measured we run a nonparametric Kendall's tau-b (Tb) 
and Spearman’s Rho correlations showing the acoustic 
phenotypes. Afterwards we assigned the morphological 
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frog sex to these resulting acoustic phenotypes separating 
sampling frogs from glacial lakes and springs and streams. 


Results 


Acoustic Phenotypes 


The canonical discriminant analysis (Figure 3A) 
showed a Wilks’ Lambda of 0.037 (Sig. 0.00). The intra 
class correlations between the discriminant variables and 
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the typified canonical discriminant functions showed the (column 1), periodogram (column 2) and spectrogram 
higher values function 1 with the FO (r=0.991) and (column 3) of each a representative pulse of each of the 
function 2 with the Number of pulses (r=0.961). We four types defined in the study (Figure 4) and their 
obtained four types of pulses defined by an oscillogram acoustical characteristics are shown in Table 1. 
Pulse type Variables 1 2 3 4 
Distortion -7.25 -5.7 -8.98 -45.9 
FO harmpow (dB) -14.66 -3.85 -19.4 -1.4 
FO harmfreq (Hz) 150.19 811.8 876.2 2011.7 
Lenght (ms) 67.61 39.63 50.89 26.37 
Number of peaks 10 20 2 73 
Peaks/ms 0.15 0.5 0.04 2.77 


Table 1: Numerical values and identifiers of the example types of pulses shown in Figure 4. Distortion (thd Matlab 
function), dB and Hz of the fundamental frequency (FO), pulse length, number of peaks and peaks per ms. 


= 
oon 


Figure 4: An oscillogram (column 1), periodogram (column 2) and spectrogram (column 3) of each a representative 
pulse of each of the four types defined in the study. 


Types of pulse 


1A A 


S 
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Among these four classes of call we found not 
significant differences between class 2 and 3 (Table 2). 
Therefore we grouped them together in one class of call, 


having a final number of three classes of call or acoustic 
phenotypes. 


Correlations N = 4 (Ratio of pulses type) 


Tau_b de Kendall 


Clases of call Class 1 Class 2 Class 3 Class 4 
see i C. of c. 1 0.667 0.667 -0.183 
Sig. bilateral 0.174 0.174 0.718 
Class 2 C. of c. 0.8 1 1,000* -0.548 
Sig. bilateral 0.2 0.01 0.279 
Class 3 C. of c. 0.8 1,000* 1 -0.548 
Sig. bilateral 0.2 0.01 0.279 
Class 4 C. of c. -0.316 -0.632 -0.632 1 
Sig. bilateral 0.684 0.368 0.368 


Rho de Spearman 


Table 2: Results of the correlation tests between each pair of the groups created according to the proportions of pulse 
types. The upper part of the table shows the results of Kendall's Tau_b, while the lower part shows the results of 
Spearman's Rho. 


When we assigned the morphological frog sex to the 
three resulting acoustic phenotypes we found that the 
frog males present two acoustical phenotypes both in the 


glacial lakes and springs and streams, whereas frog 
females presented one (Table 3). 


SVL Ratio Theoretical 
Sampling | Acoustic Frog Altitude M:F soiat ers 
sites phenotype Etog sex number |Mean (cm) (m a.s.1.) B2h (P1+P2/P Fata Won! | LeMpeLabine 
tees | | 3) (Kcal/cm2 (°C) 
year) 
P1 M 10 6.1 1636- 
: P2 M 11 6.4 1709 
Spiiigs ane 2.6:1 | 10to65 9 to 13 
streams P3 F 8 5.9 
Total 29 
P1 M 4 5.8 2213 
; P2 M 17 5.7 -2286 
Glacial lakes P3 F 5 46 25:1) 4.2:1 76 to 100 20 to 32 
Total 26 


Table 3: Breakdown of the samples taken from glacial lakes and from springs and streams, noting the number of 
individuals of each acoustic phenotype, morphologic frog sex, frog number, average and standard deviation of snout-vent 
length (SVL) in cm, sampling altitude interval, ratios of non-canonical males (P2 assimilated to “pirate” males) to 
canonical (P1 assimilated to “parental” males) and of males (P1+P2) to females (P3), theoretical solar radiation in 
Kcal/cm2 year and water temperature in °C taken during the last week of July and the first week of August from 12 to 14 
solar hours under sunny climatic conditions. 


Environmental Variables 


We compared the theoretical solar radiation and 
water temperature among our sampling sites which were 
gathered in two groups, low-and high-altitude, springs 
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and streams and glacial lakes respectively (Table 3). Both 
the theoretical solar radiation (kcal/cm2) and the water 
temperature (2C) in glacial lakes had higher values than 
those found in springs and streams (Figure 5). 


Copyright© Miramontes-Sequeiros LC, et al. 


International Journal of Zoology and Animal Biology 


42247'15-8"N 0225'45.9"W IbonEspelunciecha 


42248'49.3"N 0216'56.7"W 


Ibon de las Ranas B 


90 Kcal/cm2 year 


Figure 5: Examples of the two extreme sampling sites 
A and B (low- and high-altitude respectively) 
describing the calculation of the theoretical solar 
radiation (Kcal/cm2 year) in function of the slope 
exposition (in green), limiting sector (in yellow) and 
degrees. 


Acoustic and = Environmental 


Variables 


Phenotypes 


When comparing the acoustic phenotypes occurrence 
between both high- and low-altitude we found the 
phenotype 1 and 2 (P1 and P2_ respectively), 
corresponding to males, showed a ratio 1:1 (P2/P1) in 
low-altitude and low theoretical solar radiation and water 
temperature, whereas they had a 4.25:1 (P2/P1) ratio in 
high-altitude, high theoretical solar radiation and water 
temperature. The acoustic phenotypes male-female ratio 
(P1+P2/P3) increased with altitude, theoretical solar 
radiation and water temperature, from 2.6:1 to 4.2:1 
(Table 3). 
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Discussion 


Our results showed the existence of three acoustic 
phenotypes recognizable in the release call. The 
phenotype 1 and 2 were assigned to morphological male 
adult frogs, and phenotype 3 to morphological female 
adult frogs. In addition, the phenotype 1 (“parental” 
males) consisted of the more canonical males, which 
exhibit their vocal sacs when singing under-water, which 
emitted almost exclusively, type 1 pulses, showing a very 
tight and consistent clustering as expected from a trait 
subject to heavy selective pressure. The phenotype 2 
(“pirate” males) corresponded to individuals with male- 
like vocalizations but with different degrees of divergence 
from the canonical vocalization and closer to females, 
which did not exhibit the vocal sacs when singing under- 
water like females did. Intermediate individuals 
outnumbered both canonical males and _ canonical 
females, and they were much more abundant among 
individual sampled from high-temperature glacial lakes 
than from those taken from springs and streams. Given 
the presence of this intermediate individuals, along with 
the high male-to-female ratio observed in our study 
populations and the existence of "pirate” males that are 
unable or unwilling to attract females but fertilize 
clutches we _ propose that these high-mountain 
populations are suffering from a masculinizing effect of 
temperature that is exacerbated among those frogs that 
go through metamorphosis in high-temperature waters. 


Any normal population with exclusively genetic sex 
determination (like endotherm vertebrates) is expected 
to have close to a 1:1 male-to-female ratio, while the 
studied population has an overall 4:1 male-to-female ratio 
which is supported by deviations of sex ratios previously 
described in R. temporaria [26]. The mechanisms of sex 
determination in R. temporaria remain unclear. Their sex 
chromosomes are weakly differentiated [27,28]. And 
although some sex-linked genetic markers have been 
found Matsuba, et al. [29] others have found families and 
subpopulations in which males and females are 
genetically identical [30]. Other recent studies have found 
large differences in the extent to which sex is linked to 
genetic markers in different families and regional 
populations [31]. 


As part of the ongoing climate change, temperatures 
have been rising in the Central Pyrenees [32,33] and 
episodes of unusually high temperatures have become 
more frequent and [34]. Our study area consists of high- 
mountain valleys that are mostly free of organic 
contaminants, and the only noticeable human impact 
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comes from hikers and the occasional camper. Therefore, 
it is unlikely that any imbalance on the sex ratios is the 
consequence of environmental contaminants. Moreover, 
as a consequence the short activity period in high 
mountain and the presence of snow during spring, frog 
reproduction occurs in June and tadpole metamorphosis 
in August. Juveniles are thus exposed to the highest peak 
water temperatures of the year. Our results might suggest 
that the reduction on the proportion of females in our 
study populations is the result of rising temperatures 
having an increasingly masculinizing effect on developing 
frogs. This is reinforced by the presence of a much larger 
imbalance among those frogs that reproduce in lakes, 
which migth have been exposed to water temperatures of 
up to 32°C. This masculinization effect would result in 
males with slight anatomical and behavioral differences 
when compared to more canonical males, which would 
explain the differences in patterns and unusual behaviors 
such as male “piracy” [4]. 


Conclusion 


We differentiated three acoustic phenotypes present 
in the release call of Rana temporaria from the Central 
Pyrenees, as a _ result of studying the acoustic 
characteristics vocalization of the pulses instead of the 
entire call [35,36]. 


Two acoustic phenotypes corresponded to 
morphologic male frogs and one to morphologic female 
frogs. One of the male phenotypes (P2) might be 
represented by the “pirate” males as a result of the 
masculinization effect exacerbated by the highest water 
temperatures found in the glacial lakes. Our study 
population was acoustically homogeneous and we could 
not differentiate the occurrence of the clades mentioned 
in Teacher, et al. [2]. 
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